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ABSTRACT 

We cross correlate Gamma-Ray Bursts and X-Ray selected clusters of galaxies at 
z < 0.45. We find a positive 2(t signal for the angular cross-correlation function 
Wbc{0) on scales 6 < 3° between short GRBs and clusters. Conversely, no correlation 
is found between clusters and the population of long GRBs. The comparison with 
the cluster autocorrelation function shows that short GRBs do not trace the cluster 
distribution as not all short GRBs are found in clusters. A higher signal in Wbc{0) 
is found if we only consider the cluster population up to z=0.1. By comparing the 
short burst autocorrelation function with model predictions we then constrain short 
bursts to mostly originate within ~ 270 Mpc (i.e. z < 0.06). Our analysis also reveals 
that short GRBs are better correlated with "normal" galaxies. The double compact 
object merger model for short GRBs would associate them preferentially to early-type 
galaxies but the present statistics do not allow us to exclude that at least a fraction 
of these events might also take place in late-type galaxies, in agreement with recent 
evidences. 

Key vifords: gamma-ray: bursts — methods: data analysis 



1 INTRODUCTION 

■ The population of Gamma Ray Bursts (GRB) presents a 
bimodal duration distribution (Koveliotu et al. 1999): short 

i: GRBs (S-GRBs - lasting less than 2 s) and and long GRB 

■ (L-GRBs - lasting more than 2 s). Further support to this bi- 
, modality comes from their different spectral properties with 

short GRBs being harder than longer events (Tavani 1998; 
Ghirlanda, GhiseUini & Celotti 2004). 

One of the most accredited models associates short 
GRBs to the coalescence and final merger of two compact ob- 
jects in a binary system (Narayan, Paczynski & Piran 1992; 
Ruffert 1997). Binary systems evolve on different timescales 
(Voss & Tauris 2003) and a testable consequence of this 
model is that short GRBs should be found preferentially 
in evolved (early) type galaxies. Recently, Swift (Gehrels 
et al. 2004) and Hete-II (Lamb et al. 2004) detected for 
the first time the X-ray and optical afterglow emission of 
three short GRBs and in two cases a redshift was measured: 
GRB 050724 (Covino et al. 2005) was found to be associated 
with an elliptical galaxy at z = 0.257 (Berger et al. 2005) 
whereas GRB 050709 (Butler et al. 2005) was found to be 
associated with a blue dwarf galaxy at redshift z = 0.16 
(Covino et al. 2005, Hjorth et al. 2005) 

A third interesting case is represented by GRB 050509B 
(Hurkett et al. 2005), whose afterglow was detected only in 



X-rays, was found to be spatially coincident with a giant el- 
liptical galaxy at 2 = 0.2248 (Bloom et al. 2005) belonging 
to the cluster of galaxies ZwCl 1234.0-1-02916 at z = 0.2214 
(Pedersen et al. 2005). Another short GRB 050813 was found 
to be spatially coincident with a cluster at z = 0.7 (Gladders 
et al. 2005) . Gal- Yam et al. (2005) recently reported a third 
significant positional coincidence of the short GRB 790613 
with the rich AbeU cluster 1892 at z = 0.09. The possibil- 
ity of a correlation of GRBs with clusters of galaxies has 
been debated in the past (Kolatt & Piran 1996, Struble & 
Rood 1997) by studying the direct and statistical associa- 
tion of BATSE bursts with optically selected AbeU clusters, 
and opposite conclusions were reached (Hurley et al. 1999; 
Gorosabel & Castro-Tirado 1997). In these studies, however, 
the population of bursts was not separated in short and long 
events. Magliocchetti, Ghirlanda & Celotti (2003, hereafter 
MGC03) found evidence for anisotropy in the distribution 
of short GRBs and suggested that they might originate in 
galaxies distributed up to 2 ~ 0.5. Tanvir et al. 2005 found 
that at least a fraction between 5% and 25% of BATSE short 
GRBs might be in the very local universe (i.e. 2 < 0.025) 
preferentially in early type galaxies. It has also been pro- 
posed that a fraction (e.g. Hurley et al. 2005) of short GRBs 
might be the cosmological counterparts of the giant flares of 
Soft Gamma Ray Repeaters (SGR). 
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Figure 1. Sky distribution in Galactic coordinates of tlie sample 
of 497 BATSE short GRBs (filled circles) and of 763 Clusters (449 
REFLEX clusters - blue open circles - and 314 NORAS clusters - 
green open circles). The red filled circles represent the 283 short 
duration GRBs with position accuracy < 10° and |6| > 20°. 

In this letter we study the cross correlation between 
short GRBs detected by BATSE and X-ray selected clus- 
ters of galaxies covering the entire sky (±20° above and 
below the Galactic plane). We find that a positive corre- 
lation indeed exists between the two population of objects 
on small angular scales. We investigate if short GRBs trace 
the cluster distribution and also test if they preferentially 
correlate with early type galaxies. In the following we as- 
sume a "standard" cosmological model with Q,m = 0.3 and 
= ft = 0.7. 



2 GRBS AND CLUSTERS SAMPLES 

BATSE detected more than 2000 GRBs during its nine years 
activity (Paciesas et al. 1999). From the on-line sample of 
BATSE-GRBs^ we extracted 497 short events with dura- 
tion < 2 s and 1540 long GRBs with duration > 2 s. In 
order to correlate the GRB sample with a sample of clus- 
ters of galaxies with a sky coverage as wide as possible, we 
considered the REFLEX (Boehringer et al. 2004) and the 
NORAS (Boehringer et al. 2000) surveys. The ROSATF.SO 
Flux Limited X-ray (REFLEX) galaxy cluster survey con- 
tains 449 clusters with a maximum redshift z ~ 0.451 and it 
is flux limited to 3x10"^^ erg cm^ s"^ in the [0.1-2.4] keV 
band. It covers the Galactic latitudes |b| > 20° and decli- 
nation 5 < 2.5° and it excludes ~ 324 deg^ of sky around 
the LMC and SMC. The total area covered is 4.24 sr (i.e. 
~ 34% of the entire sky) and the catalog completeness is 
> 90% (Boehringer et al. 2004). The Northern ROSAT All- 
Sky galaxy survey (NORAS) contains 484 clusters (including 
the supplements to the original survey catalogue - Bohringer 
et al. 2000) with measured redshift up to 2 ~ 0.459. It cov- 
ers the Galactic latitudes |fe| > 20° and declination 5 > 0°. 
We are aware that the NORAS sample is only 82% com- 
plete with respect to the REFLEX survey and we tested 
our results for this difference (Sec. 3) . In combining the two 
surveys, given the higher level of completeness of REFLEX, 
we cut NORAS at 5 > 2.5°, to avoid superpositions, and to 
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Figure 2. Angular cross-correlation (red filled circles) between 
short GRBs and X — ray selected clusters of galaxies . Also shown 
is (open stars) the signal obtained with BATSE long GRBs . 
The long-dashed line represents the autocorrelation function of 
clusters, while the solid and dotted lines represent the cross corre- 
lation of clusters with early type and late type galaxies (see Sec. 
3 for details). 

the same flux limit of REFLEX. We end up with 314 (out 
of 484) NORAS clusters combined with the 449 REFLEX 
clusters (open green and blue circles in Fig.Q respectively). 

Since we do not have redshift measurements for the pop- 
ulation of short GRBs we have to rely on projected quan- 
tities. Note that for the cross correlation analysis we have 
considered only the 283 (out of 497) short GRBs (red flUed 
circles in Fig. distributed at Galactic latitudes |&| > 20° 
and with positional accuracy A < 10°. The latter selection 
limits the large uncertainties associated to the position of 
BATSE bursts (see MGC03 and Sec. 3). The same selection 
applied to the sample of 1540 long bursts gives 973 events. 
In Fig. we report the flnal samples of short GRBs (filled 
circles) and the two samples of clusters (open symbols). 



3 GRBS-CLUSTERS CROSS CORRELATION 

The angular cross correlation wtc{d) represents the frac- 
tional increase, relative to a random distribution, of the 
probability per unit solid angle of finding objects separated 
by an angle 9. We computed Wbc{0) for the 283 S-GRB (B) 
and the 762 clusters (C) by counting the burst -cluster pairs 
(BC) at different angular separations 6. This is compared 
with the pairs found at the same angular scale 9 between 
the burst sample and a random sample (BR) of ~ 10^ ob- 
jects distributed uniformly in the same sky portion covered 
by the clusters samples. 

We use the Hamilton estimator for the angular cross 
correlation function between short GRBs and clusters 
(Hamilton 1993): 



■Wtc{0) 



BC{e)RR(e) 



(1) 



http: / /cossc. gsfc.nasa.gov/docs/cgro/batse/ 



NbcNrr BR^e) 

where Nbr and Nbc represent the total number of pairs 
between the short burst sample and the random and clus- 
ter samples, respectively, and Nrr represents the number 
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of pairs in the random sample. The errors on wi,c{9) were 
derived with the bootstrap method (see e.g. Ling, Frenk & 
Barrow 1986): a set of 100 bootstrap catalogs, each the same 
size as the data catalog, were randomly extracted from the 
GRB catalog. The cross-correlation (Eq. 1) was computed 
for each bootstrap catalog and, for each 9, a set of normally 
distributed estimates of the correlation function is obtained. 
The variance around the mean represents the la uncertainty 
on wbc(Q)- In all our calculation we accounted for the BATSE 
sky exposure map (Chen & Hakkila 1997) and for the sky ex- 
posure of the REFLEX and NORAS samples, which however 
are quite uniform. As already shown in Fig. 1 of MGC03, the 
typical positional error associated with short BATSE bursts 
is of few degrees and the distribution of positional uncer- 
tainties of short bursts is nearly flat up to ~ 5° while that 
of long GRBs peaks at 1-2°. The large uncertainty in the 
positions of short GRBs might affect the correlation results. 
Following the same procedure adopted in MGC03 we tested 
how the cross-correlation signal changes by considering sub- 
samples of short GRBs with better positional determination 
(i.e. < 5°) and found results fully consistent with those pre- 
sented in Fig.|5| 

We find a positive 2a correlation signal (Fig. |5| - red 
filled circles) on small angular scales (i.e. ^ < 3°) between 
short GRBs and clusters while no correlation is found be- 
tween long GRBs and clusters (star symbols in Fig.|51l. 

Given the incompletness of the NORAS cluster sample 
(Sec. 2), we checked our results by computing the cross cor- 
relation between S-GRBs and the REFLEX and NORAS 
samples separately. Although the signal has larger uncer- 
tainties due to the smaller number of objects in the indi- 
vidual samples, we still find a positive correlation signal on 
small angular scales, entirely compatible with the one ob- 
tained by considering the two samples together. 



4 MODEL COMPARISON 

To provide an insight on the findings of the previous section 
we have to compare our data with results available in the 
literature. We use the generalization of the Limber equation 
developed by Peebles (1974) and Lilje & Efstathiou (1988), 
which relates the angular two-point cross-correlation func- 
tion Wc,i(0) to the spatial one, ^c.i(''): 



Wc,r{e) = 2 ■ 



Nc(x)dx <bi{x)x'^dx 



(2) 



where x is the comoving coordinate, u and r are related by 
the expression (which only holds in the small angle approx- 



imation, 6 « 1 rad) r 



-\- 6 X , Nc{x) is the number 



of clusters in the considered sample with distance between 
X and x + dx and the selection function for a particular class 
of extragalactic sources ^i{x) is connected to their redshift 
distribution by 



$i{x)x^dx ■ 



1 



N^{z)dz, 



(3) 



where Ws is the area covered by the survey. Note that the 
above equations have been obtained for a flat universe. 

For our analysis, we chose to consider three difi'erent 
cases for the spatial cross-correlation function: (i) = ^c,c 



(i.e. cross-correlation function coinciding with the auto- 
correlation function of clusters); (ii) ^c,i = Cc,e (cross- 
correlation between clusters and early-type galaxies); (iii) 
Cc,i = Cc,i (cross-correlation between clusters and late-type 
galaxies) . 

As for the form of the cross-correlation function to be 
plugged into equation J^J, in (i) we have used €,c,c{r) = 

with ro = I8.8//1 Mpc and 7 = 1.83, as derived from 
the analysis of the clustering properties of REFLEX clusters 
(Collins et al. 2000), while for the other two cases (ii,iii) we 
have written S,c,i{r) = ^c,c bc/bi, with b the bias factor. This 
latter expression can be easily derived as, theoretically, the 
cross-correlation function of two populations of extragalactic 
sources can be expressed as = bjbi^uM, with ^dm auto- 
correlation function of the underlying dark matter (see e.g. 
Mo, Peacock & Xia, 1993). Under the assumption of scale- 
independence, bi — as,i/as, where as is the r.m.s. of matter 
fluctuations on a scale 8/h Mpc, and ag.i can be obtained 
from measurements of the two-point autocorrelation func- 
tion for a chosen class of sources as as,i = [c-y. (''o,i/8)''*]'''^, 
with = (3_-y.)(4_^7)(6-7i) (Peebles, 1980). Values for 
the correlation length ro,i and the slope 7i for the popula- 
tions of early-type and late-type galaxies have been taken 
from studies of the clustering properties of 2dF galaxies 
(Madgwick et al., 2003) and, for as ~ 0.9 as the latest re- 
sults from CMB data seem to indicate (Spergel et al. 2003), 
we obtain be — 4.05, be — 1.54, bi = 0.97 respectively for the 
bias factor of clusters, early-type and late-type galaxies. 

Ne{x) in equation ^ has been derived from the redshift 
distribution of REFLEX clusters. Note that in this case we 
have also made the choice Ne{z) = Ni{z), i.e. we have as- 
sumed galaxies of both types to follow the redshift distribu- 
tion of REFLEX clusters. This choice was forced by the lack 
of a statistically significant spectroscopic sample of galaxies 
of different type which probes their redshift evolution from 
the local universe up to z ~ 0.45, maximum redshift of the 
clusters analysed in this work. 



5 RESULTS 

Direct comparison with the data shows that the distribu- 
tion of short-GRBs (red filled circles in Fig. |5J does not 
trace that of clusters (long dashed line in Figl^J, i.e. there 
is not a strong correspondence between clusters and short- 
GRBs, as not all short-GRBs are found to reside in clusters. 
If one instead compares the results with the cluster-galaxies 
cross-correlation function (solid and dotted lines in Fig. |21 
respectively for early- type and late- type galaxies), a much 
better agreement is obtained. 

This implies that short-GRBs exhibit similar clustering 
properties of "normal" galaxies, i.e. that they are present 
in these classes of sources. Unfortunately, the size of the 
error-bars does not allow to discern between early-type and 
late-type galaxies, preventing to associate short-GRBs with 
a particular class of galaxies. 

The local population of short bursts might be "con- 
taminated" by a more distant (i.e. z > 0.5) and isotropic 
population of objects which reduces the correlation signal 
with clusters at z < 0.5. Under this hypothesis we should 
expect the cross correlation signal to decrease when selecting 
sub-samples of clusters at lower redshifts. 
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Figure 3. Cross correlation between short GRBs and redshift 
limited samples of clusters of galaxies {z < 0.1,0.15,0.2). The 
solid (dotted) line represents the cross correlation between clus- 
ters and early type (late type) galaxies at redshifts z < 0.1. 



Figure 5. Two point angular auto-correlation function for the 
short GRB sample with positional error A < 10° . The lines rep- 
resent different autocorrelation functions for early and late type 
galaxies with limiting redshift of 0.1 and 0.05. 
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Figure 4. Cross correlation function of the sample of short GRBs 
and the clusters of galaxies sample separated in two luminosity 
sub-samples of roughly equal total number of objects. The open 
symbols represent the galaxy clusters with X-ray luminosity < 
1.4 X 10*4 erg/s (< 1.9 X lO"* erg/s) for the REFLEX (NORAS) 
sample. 

We tested for this possibility by cutting the RE- 
FLEX-I-NORAS cluster sample at different redshifts. At 
variance with what previously expected we find (Fig. 1^1 
a stronger correlation signal between S-GRBs and the 420 
clusters with 2: < 0.1 than what found with the 693 clusters 
distributed out to z = 0.2 or out to z ~ 0.4. We also tested 
for any dependence of the cross-correlation signal on the 
cluster luminosity. In order to do this, we have divided the 
REFLEX and NORAS cluster samples in two sub-samples 
by considering the median values of their X-ray luminosity 
distribution: Lx ~ 1.94 x 10*'' erg s"^ and Lx ~ 1.37 x lO""* 
erg s~^ for the NORAS and REFLEX surveys, respectively. 
This selection corresponds to have a roughly equal number 
of clusters (~424) in the two luminosity samples. As it is 



shown in Fig. 2] we find a higher signal with sub-luminous 
clsuters with respect to the signal found with more luminous 
clusters^. 

Since the cluster surveys considered in our work are 
both fiux-limited, one has that in general more luminous 
sources will mainly be more local than sub-luminous ones. 
What the data then shows is that short-GRBs might prefer- 
entially inhabit low-redshift {z ^ 0.1) systems. In fact, this 
evidence is furtherly strengthened by the comparison of our 
results with the angular cross-correlation function of clus- 
ters and galaxies, obtained as in equation (|5J by applying a 
redshift cut Zmax = 0.1. One can see (Fig. that there is 
a very good match between data and "models". Note that 
in this case the redshift distribution adopted for early-type 
and late-type galaxies has been taken in a self-consistent 
way from the 2dF Galaxy Redshift Survey (Madgwick et al. 
2002), complete at least up to 2: ~ 0.15. 

Further evidence on the proximity of short-GRBs can be 
obtained by comparing the projected auto-correlation func- 
tion for this class of sources as obtained by MGC03, with 
the autocorrelation function of both late- type and early-type 
galaxies evaluated in a fashion similar to equation ^ up to 
some maximum redshift Zmax. This is done in Fig. |3 where 
the dotted (long dashed) line represents the autocorrelation 
function of late (early) type galaxies at z < 0.1. The solid 
(dashed) line is instead that for z < 0.05. This comparison 
with the data also seems to hint to short-GRBs to be more 
local than z < 0.1, possibly being distributed only up to red- 
shifts z — 0.05. Again we note that, with the available data, 
it is not possible to discern if short-GRBs are mainly hosted 
in early-type or late-type galaxies. Better data (probably 
samples of size about three-to-five times the one considered 
in this work) are needed to draw more firm conclusions. 

^ Given the higher signal found with clusters at 2 < 0.1 we also 
tried to apply both the luminosity and redshift cuts. With these 
combined cuts our results are strongly affected by the paucity of 
clusters in the more luminous bin (e.g. only 78 with L > L, and 
z < 0.1) which does not allow us to draw any conclusion. 
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At face value, these results point to short GRBs to be lo- 
cal and preferentially associated to less dense environments 
(i.e., by extrapolation, normal galaxies) rather than to those 
identified by the clusters themselves. This last conclusion is 
reached as less X-ray-luminous clusters are generally associ- 
ated to less massive systems (see e.g. Borgani et al., 2002). 

A consequence of these results is that, given the average 
fluence (integrated > 25 keV) of 4.3x10"^ erg cm~^ s"'^ of 
our sample of short GRBs, we estimate a typical isotropic 
equivalent energy of ~ 2.4 x 10** erg, by assuming z — 
0.05. Another consequence of having short GRBs at lower 
redshift than previously thought is the increase of their rate 
as recently pointed out by Guetta & Piran 2005, Nakar et 
al. 2005 and Gal- Yam et al. 2005. 



6 CONCLUSIONS 

By cross correlating BATSE bursts and (REFLEX -I- NO- 
RAS) X-ray selected samples of clusters (distributed out to 
z ~ 0.45) we find that short duration (Tgo < 2 s) GRBs are 
correlated with clusters while we do not find any correlation 
with the population of long GRBs. By comparing the short 
GRBs-cluster correlation signal with the cluster-cluster au- 
tocorrelation signal we can exclude that short bursts trace 
the cluster distribution. Instead, through the comparison 
with the cluster-galaxy correlation functions we conclude 
that short GRBs are associated to "normal" galaxies. 

We explored the hypothesis that short GRBs are lo- 
cal events. In fact, we find a higher cross-correlation signal 
with low redshift clusters (i.e. up to z < 0.06) or with sub- 
luminous clusters. What further ly supports our conclusion 
is (i) the similarity of the short GRBs-clusters and local 
galaxies-clusters cross correlation and (ii) the short GRBs 
auto-correlation, which is similar to the auto-correlation 
function of local [z <0.1) galaxies. The present statistics 
do not allow to exclusively associate short GRBs with early- 
type galaxies as expected if they are produced - in the double 
compact merger scenario - by an old stellar population. 

These results represent a challenge because on one side 
they seem to contradict the, still very few, redshift mea- 
surements of short GRBs which place them at 2 = 0.16 
(GRB 050709 - Covino et al. 2005, Hjorth et al. 2005) and 
z = 0.25 (GRB 050724 - Berger et al. 2005) and on the other 
side they predict a typical energy of ~ 10''* erg for short 
GRBs. The apparently different (lower) redshift of BATSE 
short bursts (as found with our analysis, but see also Tanvir 
et al. 2005) with respect to the few Swift measured redshifts 
might be due to one (or to a combination) of several ef- 
fects: (i) an unknown bias of Swift to detect systematically 
larger redshifts (as also shown for the population of long 
GRBs) with respect to BATSE; (ii) a possible contamina- 
tion of the BATSE short burst population by extragalactic 
SGR (as suggested by Hurley et al. 2005), although direct 
positional searches (Nakar et al. 2005) and spectral analy- 
sis (Lazzati et al. 2005) hardly support this scenario; (iii) 
a complex luminosity function (e.g. Guetta & Piran 2005). 
For these reasons the redshift distribution of short bursts 
still represents an open issue which needs more direct red- 
shift determinations as well as a better understanding of the 
possible selection effects. 
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